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ABSTRACT 

Context. NGC 6611 and its parental cloud, the Eagle Nebula (M16), are well-studied star-forming regions, thanks to their large 
content of both OB stars and stars with disks and the observed ongoing star formation. In our previous studies of the Eagle Nebula, 
we identified 834 disk-bearing stars associated with the cloud, after detecting their excesses in NIR bands from J band to S.Oyum. 
Aims. In this paper, we study in detail the nature of a subsample of disk-bearing stars that show peculiar characteristics. They appear 
older than the other members in the V vs. V - I diagram, and/or they have one or more IRAC colors at pure photospheric values, 
despite showing NIR excesses, when optical and infrared colors are compared. 

Methods. We confirm the membership of these stars to M16 by a spectroscopic analysis. The physical properties of these stars with 
disks are studied by comparing their spectral energy distributions (SEDs) with the SEDs predicted by models of T-Tauri stars with 
disks and envelopes. 

Results. We show that the age of these stars estimated from the V vs. V - I diagram is unreliable since their V - I colors are altered 
by the light scattered by the disk into the line of sight. Only in a few cases their SEDs are compatible with models with excesses in 
V band caused by optical veiling. Candidate members with disks and photospheric IRAC colors are selected by the used NIR disk 
diagnostic, which is sensitive to moderate excesses, such as those produced by disks with low masses. In 1/3 of these cases, scattering 
of stellar flux by the disks can also be invoked. 

Conclusions. The photospheric light scattered by the disk grains into the line of sight can affect the derivation of physical parameters 
of Class II stars from photometric optical and NIR data. Besides, the disks diagnostic we defined are useful for selecting stars with 
disks, even those with moderate excesses or whose optical colors are altered by veiling or photospheric scattered light. 
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1. Introduction 

The evolution and the physical properties of circumstellar disks 
around young pre-main sequence (PMS) stars is currently a hot 
topic of astrophysics. These objects, in fact, provide the environ- 
ment where planetary formation takes place, and we are far from 
fully understanding the mechanisms and conditions that make 
the disks evolve toward planetary systems. 

With current telescopes, circumstellar disks can be spatially 
resolved only in very nearby star-forming regions (less distant 
than a few hundred parsecs). In more distant targets, photometric 
observations allow selecting a large number of PMS stars with 
disks by detecting their infrared excesses, with respect to the 
photospheric emission, from the disk itself. In these cases, the 
physical parameters of disks can be studied through the analysis 
of their infrared emission using models of stars with disks, and 
possibly outer envelope. In particular, the observed spectral en- 
ergy distributions (SEDs), i.e. the fluxes emitted by the sources 
in various photometric bands vs. their central wavelengths, can 
be compared to models predicting SEDs, in order to estimate the 
physical properties of the systems. This is an indirect method 
of studying the parameters of stars with disks, strongly depend- 
ing on the reliability of the chosen models. However, in some 
cases the estimates of the physical parameters of accreting disk- 
bearing stars obtained with SED analysis has been successfully 
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compare d with the values deriv ed with other independent meth- 
ods (e.g. lRobitaille et al.ll2007l) . 

Since the public ation of the ini t ial mo dels of astrophysica l 
accretion disks (i.e. IShakura et all Il978l and iPringlel 1 198 lb . 
several authors developed models of circumstellar disks around 
PMS stars. Recently, a deta iled model of the structure of thes e 
disks has been developed by D'Alessi o et al](ll998lll999ll2001h . 
A full summary of the physical characteristics and the simplify- 
ing assumptions of this model can be found in D' Alessio et al.l 
(i2005i) . The available model grid comprises 18 different central 
stars with four disk models each (characterized by four different 
accretion rates). 

More recentlv. fRobitaille et alj(l2006l) have presented a set of 
models that covers an extensive parameters space: 20 000 young 
stellar objects (YSO) models, each at 10 different inclination an- 
gles, for a total of 200000 distinct SEDs. These models take the 
emission of the collapsing envelope into account and allow sep- 
arating the photospheric light directly observed from the light 
scattered by the disk into the line of sight. This feature of these 
models is crucial for this work. 

In this paper we analyze the SEDs of a subsample of can- 
didate disk-bearing stars associated to the Eagle Nebula (M16). 
This star-forming region, and the associated clu ster NGC 661 1, 
has be en the subject of our precedent studi es: Guarcello et al.l 
(2007), hereafter GPM07, Guarcello et al.' ('2009"), hereafter 
GDM09, and Guai'cello et al. (2010), hereafter GMP 10. This pa- 
per is organized as follows. In Sect. |2] we summarize the param- 
eters of NGC 661 1 and M16 found in our preceding studies; in 
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Sect. [3]we discuss the disk diagnostics used in our previous pa- 
pers; in Sect. |4]we analyze and discuss the SEDs of candidate 
members with disks that in the V vs. V - 1 diagram appear older 
than the other cluster members; in Sect. |5] we study the SEDs of 
the candidate members with disks with one or more IRAC colors 
at photospheric values; the final results are summarized in|6] 

2. Cluster parameters and analyzed data 

In our previous studies of the young open cluster NGC 66 11 and 
its parental cloud Ml 6, we compiled a multiband catalog of this 
region including 190684 sources, using 

- optical data in BVI bands obtained with the Wield Field 
Camera (WFI) at ESO/2.2m fro m observations which are 
part of the ESO Imaging Survey (iMomanv et al.L 1200 Ih and 
which have been reduced and analyzed in GPM07; 

- NIR data in JHK bands obtained fro m the 2 Micron All Sky 
Survey (2MASS, ICutri et al.L l2003h Point Source Catalog, 
analyzed in GPM07; 

- NIR data in JHK ba nds taken from the Galactic Plane Survey 
dLucas et al.L |2008') of the United Kingdom Infrared Deep 
Sky Survey (UKIDSS) analyzed in GMPIO; 

- Spitzer/IRAC data in four bands centered on 3.6//m, 4.5/im, 
5.8//m and 8.0/zm obtained from the Galactic Legacy 
Infrared Mid-jPlane Survey Extraordinaire (GLIMPSE; 
iBeniamin et"al]|2003h . analyzed in GDM09; 

- X-ray data obtained from three Chandra/ACIS-I observa- 
tions: one with 78 Ksec of exposure time (P.I. Linsky) 
centered on NGC 6611 and two of 80 Ksec, analyzed in 
GMPIO (P.I. Guarcello), centered on the pillar of gas called 
"Column V" (iOliveirai l2008i) and on a dense region of the 
molecular cloud. 

The WFI, 2MASS, and UKIDSS data are available in a field 
of 33' X 34' centered on NGC 6611 (15.8 x 16.2 parsec at the 
distance of the cloud); also the IRAC data are available in this 
field, with the exception of a region northwest of about 7' x 10'. 
X-ray data cover about 3/4 of this studied field. Data at even 
long er waveleng t hs of M16 are provided by the MIPSGAL sur- 
vey dCarev et al. ,l2009b. but th e observation is affected by an in- 
tense and patchy background (Flagev et al., 2009) in correspon- 
dence with NGC661 1, which complicates the estimation of point 
sources flux. For these reasons, these data are not studied in the 
present work. 

In our previous works we evaluated several physical prop- 
erties of M16 and its stellar population, such as distance, age 
spread, and average extinction. In these works, the influence of 
the OB members of NGC 66 1 1 on the evolution of circumstellar 
disks and the star formation in the nebula was also studied. Table 
[U summarizes some of the physical pro perties of NGC 66 1 1 and 
M16 that we found. See lOliveiral (l2008l) for a complete review 
of the Eagle Nebula. 

In our previous studies we identified 834 candidate disk- 
bearing stars. Among them, there are several stars with inter- 
esting color peculiarities (discussed in the following sections). 
In this paper we s tudy their properties u sing the SED fitting 
tooQ developed bv iRobitaille et al.] (l2006l) . In our analysis, for 
each studied star, we select the models that best reproduce 
the observed SEDs, using the compatibility criterion defined in 
[Robitaille et al. 2007 (i.e. those for which the reduced chi-square 
satisfies the relation ~Xbest - 'where;^'^ 
fit model). 



Table 1. Physical properties of NGC 661 1 and M16 obtained in 
GPM07, GDM08, and GMPIO. 



tf^^^^ is that of the best- 



Parameter 



Distance 
Median age of NGC 6611 
Age range of stars in Ml 6 
Average DF* in NGC 6611 
Average DF* far away from OB stars 
Average DF* close to OB stars 
Average Ay in NGC 661 1 
Rv 

Core relaxation time 



Value 
1750 parsec 
~ 1 Myear 
-3-0.3 Myears 
36% ± 1% 
42% ±3% 
28% ± 3% 
2.6'" 
3.3 
5.2 Myears 



*DF stands for Disk Frequency 

3. Reddening-free color indices used for disk 
selection. 

To select stars with infrared excesses, in our previous studies 
we used a set of color indices {Q indices), which are defined to 
be independent of reddening ( fPamiani et all 120061) . The set is 
composed of nine Q reddening-free color indices, defined with 
the bands used in our previous studies. Eight indices compare 
the infrared colors with V - I, assumed to be representative of 
photospheric emission: 



QviAB = (V - /) - (A - B) x Ev-i/Ea-b 

and one compares J - H with H - K, i.e. 
QjHHK ^(J-H)-(H-K)x Ej^hIEh-K 



(1) 



(2) 



available at |http://caravan.astro.wisc.edu/protostars/index.php| 



In these equations, A - B is an infrared color (/ - J, 
J - H, J - K, H - K, J - [3.6], J - [4.5], J - [5.8] and 
J - [8.0]); Ev-i, Ej-H, Eh-k, and Ea-b are reddening in the 
corresponding colors, which have been ca lc ulated f r om th e red- 
dening laws o f Rieke & Le bofskvl (Il985h . iMafliisI (Il990h . and 
Munari & Carrara (i 1996 ). Hereafter, we refer to the whole list of 
indices with Q or Qabcd, and call Qimass and Qukidss the in- 
dices defined with VIJHK photometry and IR data respectively 
from 2MASS/PSC and UKIDSS/GPS, while Q,rac and Qvij[sp] 
are those involving IRAC bands. 

Figure [1] illustrates the case of the Qvu[3.6] vs. J - [3.6] di- 
agram used to select stars with excesses in [3.6] in GDM09. In 
this figure are plotted both the color and index of all the stars 
falling in the surveyed field with crv-i and cry_[3 gj smaller than 
0.15"'. The reddening vector is horizontal since the Q index is in- 
dependent of extinction. For this reason, in the Qabcd diagrams, 
like that shown in Fig. [1] reddening shifts the star's positions 
rightward along the x direction, from the locus of "normal stars" 
to the locus of non-classifiable stars (reddened stars and/or stars 
with excesses). In Fig.[T] if 7 - [3.6] is redder than the photo- 
spheric values since there is an excess in [3.6], the index de- 
creases while J - [3.6] increases, shifting stars downward to the 
locus of stars with excesses. In this way, stars with excesses and 
reddened photospheres can often be separated, even if a num- 
ber of stars remain unclassified since they fall in the locus where 
we cannot distinguish between reddening or excesses. The use 
of 9 2 indices to detect excesses in every single infrared band 
reduces the number of unclassifiable stars. In the context of this 
paper, it is important to note that the indices decrease even if 
V -I decreases. This can also select stars having V -I bluer than 
the photospheric values (such as when there is intense optical 
veiling). 
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Fig. 1. Qvij[3.6] vs. 7 - [3.6] diagram of the stars in the studied 
field (points). The loci of normal stars, stars with excesses in 
[3.6] and that of stars for which the effects of reddening and 
excesses cannot be discerned are shown. Stars with excesses in 
[3.6] are marked by circles. The reddening vector is also shown. 



For each Qabcd index, the stars with excess in D band are 
defined as those with the index significantly (i.e. by more than 
3 times the error in the Q index) smaller than the lower limit of 
photospheric Q values (the horizontal dotted-dashed line in Fig. 
[T]l. These limits have been defined as the lower boundaries of 
the loci, in these diagrams, of the stars with normal colors. For 
QiMASs and Qukidss , the locus of normal stars has been defined 
by using th e photospheric colors obtained from the isochrones of 
ISiess etaLldioOOl) . with age ranging from 0.25 Myears to 10^ yrs 
and the masses ranging from 0.2-7 Mq in order to consider all 
possib le observed disk-less stars. The isochrones of ISiess et al.l 
(l2000h do not provide IRAC colors, so for Qirac the locus of 
normal stars have been defined using the sources with photo- 
spheric IRAC colors and optical emission. GDM09 have shown 
that the completeness of the selection of candidate stars with 
disks made with Q indices depends on both the mass and age of 
the disk-bearing stars. 



4. Blue stars with excesses 

The left panel of Fig. |2] shows the y vs. V - I diagram for 
stars in the surveyed region with errors smaller than 0.1™ and 
0.15™ in V and V - I respectively. Since X-rays emission of 
PMS stars is much more intense than for stars in main sequence 
(MS) phase, in this diagram the positions of stars detected in 
X-rays define the locus populated by the young stars associated 
with Ml 6, which is de limited by the 0.1-3 Myears isochrones 
(from lSiess et al.Lr2000h . In the V vs. V -I diagram, however, 90 
stars with NIR excesses, detected by the diagnostics described 
in Sect. [3] have the V - 1 bluer than the PMS locus, i.e. with po- 
sitions typical of field main sequence stars older than the other 
cluster members. Hereafter we call these stars "blue stars with 
excesses", or with the acronym BWE stars. The right panel of 
Fig. |2] shows the spatial distribution of all the stars in our multi- 
band catalog, disk-bearing stars and the stars with disks which 
appear older that the cluster locus in the V vs. V - / diagram, 
falling within the 33' x 34' region centered on NGC661 1. 



Before discussing the characteristics of BWE stars in detail, 
it is important to stress that GMP07, GDM09, and GMPIO have 
been conservative in defining the stars with excesses, in order to 
avoid the contamination of the members list as much as possible. 
We also stress that the following characteristics allow us to be 
confident about the member status of most BWE stars: 

- In the IRAC color-color diagram, the only 4 BWE stars with 
emission in all IRAC bands exhibit colors compatible with 
ClassII YSOs. 

- In the NIR color-magnitude and color-color diagrams all 
BWE stars are consistent with being cluster members, with 
the exception of 9 stars that have NIR colors typical of fore- 
ground sources and that are be excluded from the following 
analysis (since they are likely field stars). 

- In 45 BWE stars with small JHK errors, the excesses have 
only been detected with the Q indices defined with opti- 
cal and infrared photometry. This does not exclude these 
stars being mismatches between optical and NIR catalog. 
However, in 30 cases also the Qjhhk index shows a weak 
/T-band excess, although it is not significant (at 3cr), inde- 
pendently of optical photometry. In fact, the mean value of 
the Qjhhk index for the entire sample of BWE stars is equal 
to -0.23, much lower than the mean Qjhhk of all stars in the 
studied FOV (equal to 0.30), pointing toward the member- 
ship of the majority of this sample of stars. 

- Two thirds of BWE stars have infrared excesses in JHK both 
using the UKIDSS/GPS and the 2MASS/PSC photometry. 

Another relevant hint about the nature of BWE stars as 
disk-bearing YSOs arises from the preliminary analysis of 
optical spectroscopic observations. We obtained observations 
(with a total exposure time of 6 hours) of the 20 bright- 
est BWE stars with the intermediate resolution spectrome- 
ters GIRAFFE/FLAMES @VLT (Obs. ID: 083.C-0837; PL: 
Guarcello). These observations have been performed with the 
HR15N{R^11QQ0) setup, which covers the 5200-7600A spectral 
range, comprising both the Ha 6562. 8A and Lithium 6707. 8A 
lines. Among the BWE stars observed with FLAMES, 6 have 
an evident Li absorption line and 13 broad and asymmetric Ha 
line, thus showing spectroscopic evidence of their PMS nature 
(4 among these stars have both features). Among the BWE stars 
with accretion signatures, eight have a blueshifted absorption 
feature (with respect to the center of the line), which is a sig- 
nature of occurring inflow of gas onto the stellar surface; two 
stars have an asymmetric profile of the line; three stars have both 
the absorption f eature and the asymmetric profile. Following 
'Kur osawa et al.l (|2P06), these line profiles indicate both gas ac- 
cretion onto the stars and wind from the disk surface with differ- 
ent values of accretion rate, mass-loss rate, and disks inclination. 
A full account of the spectroscopic data will be the subject of a 
forthcoming paper. In the present context, they show that about 
3/4 of the BWE stars are indeed young PMS stars. 

The spatial distribution of BWE stars (right panel of Fig. |2]l 
may suggest that it is uncorrelated with the cluster. However, in- 
side NGC 661 1 BWE stars are more frequent at large distances 
from OB members, which are concentrated in the center of the 
cluster. It is possible to suppose that the physical mechanisms 
responsible for the anomalous optical colors of BWE stars are 
more prominent at large distances from the center of the cluster. 
This hypothesis is consistent with the smaller disk fraction close 
to massive stars observed in our previous studies, and it can be 
due to a lower incident UV radiation or to the dynamical evo- 
lution of the cluster. Moreover, our studies of the outer regions 
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Fig. 2. Left panel: V vs. V -I diagram of all stars in the studied FOV (gray points), X-ray sources (crosses), and BWE stars (boxes). 
The thick solid line is the ZAMS (from Si ess et al. .2000) at the distance of 1750 pc and with Ay = 1.45"'. The dashed lines are 
the isochrones at 0.1, 0.25, 1, 2.5, 3, and 5 Myrs, with the average ex tinction appropriate for cluster members (Ay - 2.6"'). The 
extinction vector is obtained from the law of 'Mun ari & Carrara (Il996h . Right panel: spatial distribution of stars in the multiband 
catalog of M16 (gray points), candidate disk-bearing stars (black dots), and BWE stars (boxes). North is up, east is on the left. 



of the Eagle Nebula revealed a large number of candidate young 
stars associated with the outer region of M16. Most of the BWE 
stars, then, could be part of this widespread population of M16. 

4.1. Hypotheses on the nature of BWE stars 

One possible explanation of the nature of BWE stars is that 
they are strong accretors, characterized by an intense veiling 
that alters their optical colors. The possibility that the position 
of accreting T-Tauri stars in the H-R diagram can be affected 
by veiling has been discuss ed by several authors. For example, 
iHartmann & KenvonI (Il990h studied 37 K7-M1 T-Tauri stars in 
the Taurus- Auriga complex, and they found that strong accretion 
(M > lO^^Moyr"') affects the positions in the H-R diagrams of 
young stars, from which it is impossible to have an accurate esti- 
mate of their age. Stars with properties similar to our BWE stars 
have been observed in other star-forming regions. For example, 
Hillenbrand (1997) obtained the V vs. V - I diagram, down to 
V ~ 22"', of the Orion Nebula Cluster, where a few tens of K-M 
cluster members are even bluer than the ZAMS at the distance 
of the nebula. This author hypothesized that these stars can be 
heavily veiled or their photometry might be contaminated by the 
nebula emission. 

Another possible explanation for the observed colors of the 
BWE stars is that they are cluster members surrounded by a 
dusty reflection nebula. In this case, the extinction at / and in- 
frared bands is different from that at V, since the optical light 
is reflected by the nebula and the optical colors are very dif- 
ferent from what is expected fr om a reddened photo sphere. 
This hypothesis was proposed by iDamiani et a n (120061) for the 
young open cluster NGC 6530 to explain stars with excess de- 
tected with Qviij, which, in the optical color-magnitude dia- 
gram, are apparently older than the X-ray emitting members of 
the cluster (i.e. like our B WE stars). To support this hypothe- 
sis, Damiani et al.l (1200 6*) pointed out that two of the stars in 
the Taurus-Auriga complex, with excess detected with Qviu, are 
Class I sources surrounded by an envelope, with a faint optical 
counterpart. This hypothesis is also supported by the direct high 



spatial resolution observations of the Hubble Space Telescope 
of stars surrounded by a refl ection nebula in the Taurus-Auriga 
cluster dPadgettet al.Lil999l) . 

It is also possible that a significant fraction of stellar ra- 
diation is scattered into the line of sight by the dust grains in 
the disk atmosphere. If the light emitted by the star is scattered 
away from the line of sight, such as while passing through the 
isotropic interstellar medium (ISM), the source is affected by 
normal interstellar extinction. Stars with disks, however, are not 
isotropic systems. The amount of light scattered into the line of 
sight by the circumstellar material can be greater than the ex- 
tincted amount of radiation. Since small dust grains, with sizes 
< lju m, scatter opti cal light more efficiently than infrared radi- 
ation (iThroop et al.L 1200 ll and references therein), if the disk 
surface is mostly populated by small dusty grains, the optical 
radiation at short wavelengths is more efficiently scattered. If a 
large amount of such radiation is scattered into the line of sight, 
the observed flux emitted by the reddened photosphere at shorter 
wavelengths increases, giving a blue V - I color The efficiency 
of this mechanism depends on the distribution of the grain sizes 
in the disk atmosphere (which changes during disk evolution) 
and on other disk properties such as the disk incUnation angle 
and its scale height. 

We cannot a priori reject the hypothesis that BWE stars are 
older than the other PMS members. There are known examples 
of classical T-Tauri stars much older than the disk survival times, 
which are usually observed or predict ed by the YSOs evolution- 
ary scenario. For example, P alla et al] (1200 5. 2007) studied some 
lithium-depleted members of the Orion Nebula Cluster, among 
which they found a Classical T-Tauri star with large NIR ex- 
cesses despite its age (estimated from Li abundance) equal to 
~ 20 Myrs, and lArgiroffi et all ( l2007h studied MP Muscae, a T- 
Tauri star with disk, low accretion (~ lO""M0/yr) and an aged 
~ 17 Myrs. 

Finally, it is also possible that BWE stars come from to mis- 
matches between optical and infrared catalogs. A random coin- 
cidence between a foreground MS optical star and an infrared 
source, in fact, can produce a source with negative Q indices 
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and "blue" optical colors. This is indeed a reasonable hypothe- 
sis for the 9 BWE stars with H - K compatible with field stars. 
However, the mean value of Qjhhk for all BWE stars (-0.23) 
implies that this sample is not dominated by mismatches. 

4.2. SED analysis of BWE stars: accretion 

The properties of the 81 BWE stars (all the 90 BWE stars mi- 
nus those with H - K color compatible with field stars) have 
been studied by SED analysis, as described in Sect. [1] The ex- 
plored parameters space is described in Robitaille et al. (2006). 
Among the analyzed BWE stars, 43 are compatible with at least 
one SED model with a good;^'^ value. Mostly thanks to the small 
photometric errors in UKIDSS bands (GMPIO), these stars are 
typically compatible with just a few models, and the member- 
ship of 9 among them is confirmed by our spectroscopic data. 
For the remaining stars, their membership status is supported by 
the properties described in Sect. |4] but we cannot discern their 
nature by SED analysis. 

As explained previously, high accretion rates, and the result- 
ing intense optical veiling, could account for the blue optical 
colors observed in the BWE stars. However, only 5 BWE stars 
have the SEDs compatible with at least one model that shows, in 
the optical bands, non photospheric excesses due to gas accre- 
tion (hereafter, we call these stars ''E-stars"). 

The analysis of thek SEDs suggests that M is on the average 
higher in the E-stars sample from all the other candidate mem- 
bers with a disk, nevertheless the observability of the excess in 
optical bands requires not only high accretion rates, but also low 
stellar temperatures. In fact, in stars with low effective tempera- 
ture, the contrast between the photospheric emission and a possi- 
ble excess in V band is higher than in stars with higher effective 
temperature. Since the SEDs of all but one E-star are also com- 
patible with models without these optical excesses, we conclude 
that the optical veiling is the correct explanation for the optical 
colors of BWE stars only in few cases. 

4.3. SED analysis of BWE stars: scattering and obscuration 
by the disk 

The optical colors of BWE stars might be alternatively explained 
with the interaction between stellar radiation and the circum- 
stellar disk. The models of iRobitaille et al.i (i2006l) allow us to 
separate the stellar radiation that is directly observed from that 
scattered into the line of sight by disk's dust grains. The amount 
of light scattered into the line of sight can significantly alter the 
stellar SED (which should be a reddened photosphere at the ef- 
fective stellar temperature), by increasing the observed optical 
flux mostly at short wavelengths. Moreover, disks seen at large 
inclinations can partially occult the stellar photosphere, so part 
of stellar radiation is totally absorbed by the circumstellar ma- 
terial. As a consequence, a reduced photospheric flux with un- 
changed colors (since they are related to the stellar temperature) 
is observed in these cases. 

Among the 43 BWE stars with good SED fitting, it predicts a 
large fraction of optical emission composed by scattered light in 
26 stars, and in 16 cases a partial occultation of the stellar pho- 
tosphere by the disk without scattering effects. In the following 
we discuss the results obtained from the analysis of the former 
26 BWE stars. 

As an example, the left panel in Fig. [3] shows the SED 
model compatible with the observed SED of 661 1-2775(0, while 

^ This ID identifies the stars in the catalog published in GMPIO 



the right panel shows the individual components of the model. 
Comparing, in the right panel, the optical scattered flux with the 
directly observed photospheric flux, it is evident that the light 
scattered by the disk into the line of sight dominates the SED 
at optical wavelengths. Therefore, the optical SED is not sim- 
ply consistent with a reddened photosphere (which is also di- 
rectly observed). The scattering increases the optical flux signif- 
icantly (with respect to the expected photospheric values) mostly 
at short wavelengths, accounting for the optical blue colors ob- 
served for this star 

In the right panel of Fig. [3] the optical part of the SED 
is dominated by scattered light. If we calculate, in a particular 
band, the fraction of observed flux that is scattered into the line 
of sight on to the total emitted flux as 

Fracscat = — (3) 

^ dir ^ scat 

we find that for almost all the BWE stars with significant, 
scattered optical light, more than 90% of the observed optical 
flux is scattered by cuxumstellar material: without scattered light 
these stars could n ot have been detected in optical bands. 

The models of IRobitaille et al.l (l2006l) do not provide infor- 
mation on the grain population in the disks (for example the av- 
erage grain size), so we cannot test whether the disks in BWE 
stars with a large amount of scattered light are mostly populated 
by small grains, as expected. However, collisional aggregation 
and settling in the disk midplane is expected to occur during the 
evolution of the disks. Thanks to these phenomena, the num- 
ber of small grains in the upper layers of the disk decreases 
during the evolution of the disk itself, thus also decreasing the 
amount of scattered optical light. The typical timescales for set- 
tling and aggregation are not well determined, and some mod- 
els (as that developed by Dullemond & Dominik, 2004) predict 
short timescales of about 1 Myear for settling to occur. However, 
this disagrees with some older Herbig stars with disk having 
been observed with a significant amount of scattered light. We 
estimated the age of BWE stars with a significant amount of scat- 
tered flux using the SED models. The results are summarized 
in Table |2] showing that the SED analysis for almost all these 
stars predicts an age under 1 Myear, younger than the mem- 
bers median age (~ 1 Myear, GMPIO). This suggests that, on 
average, a large amount of scattered light is present among the 
youngest cluster members, probably associated to a large frac- 
tion of small grains in the disk upper layers, still present at 
young ages. However, we caution the reader not to overinter- 
pret our finding since it is very difficult to assess whether this is 
a r eal result, or if it is produced by some assumption in the mod- 
els. |Dullemond&|£omiiiik( 2004) list some physical effects, not 
considered by the available disk models, which can explain the 
presence of grains with typical sizes of those in the ISM in the 
upper layers of disks older than 1 Myear. We focus our attention 
just on the large difference between the age of the BWE stars 
predicted by the models (less than 1 Myear and compatible with 
the cluster age) and the age that can be deduced from the V vs. 
V -I diagram (several tens of million of years), likely caused by 
the scattered optical light. 

As explained, among the 42 BWE stars for which it was pos- 
sible to study their SED, in 16 cases the central star is largely oc- 
culted by the disk, and there are no evident scattering effects. In 
these cases, these stars are much fainter in optical band than ex- 
pected from the stellar properties, while the V - / is not affected 
by the occultation having the values that should to be expected 
from the stellar effective temperature. SED analysis predicts that 
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Fig. 3. Left panel: the SED model (solid lines) that reproduces the observed SED (black points) of a BWE star (66 1 1-27750). The 
dashed line is the unreddened photospheric flux expected if the star should to be without disk. The right panel shows the different 
components of the SED model, marked with different kind of lines, as explained in the top of the panel. Black points mark the 
observed fluxes, which are corrected for interstellar extinction using the Ay predicted by the SED model (2.3™). Both models are 
normalized at the distance of 1 Kpc. 



Table 2. Age of the 26 BWE stars with significant amount of 
scattered flux predicted by SED fitting. 



Age(years) 


^stars 


< 10' 


2 


10^ -=- 10* 


21 


10* 10' 


3 



these stars are on average older that those with a large amount 
of optical light scattered into the line of sight, as shown in Table 
[3] where we report the age of these 16 stars predicted by the 
models. 

Table 3. Age of the 16 BWE with an obscured photosphere pre- 
dicted by SED fitting. 



Age(years) 




< 10' 


1 


10^ 10*^ 


2 


10« ^ W 


13 



The SEDs of all the analyzed 42 BWE stars (both with an 
obscured photosphere and scattering effects), but four, are mod- 
eled with highly inclined disks (with an inclination greater than 
80° with respect to the line of sight). This indicates, as expected, 
that scattering effects are more important for disks observed at 
high inclination. To verify this, we produced the optical color- 
magnitude diagrams using all the ClassII models in the database 
of Robitaille et al. (2006). Figure|4]shows the V vs. V-Zdiagram 
obtained with the magnitudes of ClassII objects predicted by the 
models. All the available models of ClassII objects are plotted 
in the upper panel. They populate mostly the PMS locus, to the 
right of the 10 Myrs isochrone, but a number of them lie also in 
the locus typical of older stars, which is leftward of the 10 Myrs 
isochrone. Only the models for ClassII objects seen at low disk 
inclination angle (6 < 80°) are plotted in the lower panel. In this 



diagram, points lie only to the right of the 10 Myears isochrone. 
The comparison between the two panels suggests that the effects 
related to a disk observed at high inclination (both scattering and 
obscuration) shift the position of the stars in the V vs. V - / dia- 
gram leftward of the locus of PMS stars. 






-1 1 2 rj 

V-] 



5 




^^ 1 

-1 1 2 3 

V-] 

Fig. 4. y vs. V - I diagrams obtained from the magnitudes o f 
ClassII YSOs predicted by the models of Robitaifle et al.'(l 2006l) . 
The line is the 10 Myrs isochrone of Siess et al. (2000). In the 
upper panel all ClassII models are plotted in the lower only those 
with disk inclination smaller than 80°. 



6 



M. G. Guarcello et al.: Disk-bearing stars observed in scattered light 



5. Candidate stars with NIR excesses and 
photospheric IRAC colors 
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Fig. 5. [3.6] vs [3.6]-[4.5] diagram of the stars in the surveyed 
region (small points). Boxes mark stars with NIR excesses se- 
lected with Q indices. The reddening vector with Ay^l S™ was 
obtained from the reddening law of Meaeath et al. (2004), while 
that with Ak-1.5'" from that of .Flahertv et al.. (2007.) 



Figure |5] shows the [3.6] vs. [3.6]-[4.5] diagram for all the 
stars in our catalog with o"p g] < 0.1™ and cr[3 6]_[4 5] < 0.15"'. 
Most candidate stars with disks lie in the region of this diagram 
with [3.6] - [4.5] > 0.2'", showing a red color, as expected. In 
Figure|5] 66 stars with excesses have [3.6] - [4.5] < 0.2'", thus 
showing a [3.6] - [4.5] color at pure photospheric value even if 
infrared excesses are detected by the use of the Qirac indices. In 
the following, we describe why a star with disk could have such 
[3.6]-[4.5] colors, while having excesses in Qirac indices, and 
then we discern the nature of these stars studying their SEDs. 

It is possible that the disk has an evacuated inner region, re- 
sulting in a weak emission at the NIR bands. Also a high disk 
inclination can self-mask the NIR emission from the inner part 
of the disk. However, such disks would be selected with Qvij[s.o], 
which can be affected by a prominent silicate band in emission. 
Another possibility is that the disk has a low mass, producing 
moderate NIR excesses with respect to the photospheric emis- 
sion. In this case the [3.6]-[4.5] color is not strongly affected by 
the emission of disks, but 2v/j[3.6] and Qvij[4.5] indices can also 
detect these moderate excesses since they compare J - [3.6] or 
7 - [4.5] to y - /. This depends on the spectral type of the central 
star, as shown by GDM09. 

Among these 66 stars, only in 33 cases the physical disk pa- 
rameters are determined well with SED analysis. The models 
compatible with the SEDs of 22 among these stars range over 
all the hypotheses listed above. A fraction of these stars (10) 
have SEDs compatible with disks having masses < IO^^Mq, ac- 
counting for moderate excesses with respect to the photospheric 
emission at [3.6] and [4.5]. It is remarkable that Q indices, since 
directly comparing NIR colors with V-I, can be effectively used 
to select these disks. The SEDs of the remaining 12 stars predict 




Stellyr Mass (M^) 



Fig. 6. Diagrams of Qjhhk evaluated from the models of 
iRobitaille et aP (l2006l) of ClassII stars with thin disk, vs. the 
mass of the central star (upper panel) and the disk mass (lower 
panel). The dotted line is the limit of photospheric emission in 
the Qjhhk index. 



an almost evacuated inner disk with a prominent silicate emis- 
sion band. 

To verify the hypothesis that the adopted Q indices can select 
objects with low-mass disks, we used the NIR photometry, in 
the 2MASS system, predicted by the models of Robitaille et alj 
(2006), to compute the expected Qjhhk index of disk-bearing 
objects. We evaluated this index for models of ClassII YSOs 
with thin and evolved disk, i.e. for which Mdisk/Ms,a,- < 10"^. 
Fig. |6] shows the diagrams of their Qjhhk vs. the mass of the 
central star (upper panel) and the mass of the disk (lower panel). 
In both panels several models have Qjhhk below the limit of 
photospheric emission, and then they would be selected with this 
index. The lower panel shows that most of these stars have disks 
with very low masses (from 10"^ to 10"^ Mq). For even lower 
masses the disk emission becomes too faint to give a negative 
Qjhhk- In the upper panel, it is shown that negative indices are 
obtained for stars with masses over ~ 3Mq, likely because less 
massive stars cannot heat the material in the disk sufficiently, 
so its emission remains weak. These diagrams confirm the effec- 
tiveness of Q indices in detecting thin disks around intermediate- 
and high-mass stars. 

The SEDs of 1 1 stars cannot be reconciled with any of the 
considered hypotheses. For these stars, the SED analysis sug- 
gests that the excesses detected with various Q indices are due 
to an alteration of the optical SED by scattered light, as for the 
BWE stars. In the following, we use the example of the source 
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6611-20515 (see Fig.|7]l to illustrate this effect. As in the case 
of 661 1-27750 (Fig.O, the optical flux of this star is dominated 
by scattered light and not by the directly observed photospheric 
emission. 

Table 4. Stellar and disk parameters of the source 6611-20515 
obtained with the SED analysis. 



Parameter 



Value 



Age (Myrs) 

M* (Mo) 
Mjisk (Mo) 

R,„ (AU) 
Inclination (degree) 
M (Mo/jr) 



< 0.2 
< 1 

1 X 10"^ ^ 1 X 10-2 
> 10 
30^60 
10-^ ^ 10-' 



The parameters of 6611-20515 obtained with SED analy- 
sis are summarized in Table |4] The excesses of 6611-20515 
are detected using Qimass^ Qukidss and, despite the normal 
[3.6] - [4.5] color (equal to 0.08), 2v/y[3.6] and Qvu^.s] indices. 
The left panel of Fig.|7]shows the observed SED of this star, with 
the 5 compatible model SEDs overplotted, while the right panel 
shows the components (i.e. the emission from star and disk) of 
the best compatible model. This model predicts a large inner disk 
radius equal to 26 AU, with the disk emission becoming signif- 
icant only for A > IQfiin, not accessible with our data, and that 
of the envelope at even longer wavelengths. The flux at shorter 
wavelengths (comprising all the bands in which the excesses 
have been detected) is emitted only by the photosphere, accord- 
ing to the observed value of [3.6] - [4.5]. However, in the right 
panel of Fig.|2l it is evident that the observed stellar flux is com- 
posed of both direct and scattered fluxes, as for the BWE stars. 
In particular, the optical/NlR emission is dominated by scattered 
light (more than 90%) for A < l/um and by direct emission at 
longer wavelengths. The boundary between the regions of the 
spectrum dominated respectively by direct and scattered flux is 
at about l/um, between WFl 1 band (O.dijjm) and J band (1 .2/im). 
As a consequence, the non photospheric Qvij[3.6] and Qvij[4.5] 
indices are produced by the different paths of the photospheric 
flux at these wavelengths, since J - [3.6] is at a photospheric 
value, while V - / is bluer due to the scattered light. The ex- 
cesses detected in this star, then, are not related to the emission 
from the inner disk (that is evacuated) but they are due to the 
scatter of a significant amount of stellar light into the line of 
sight by disk and envelope. A further confirmation of the emis- 
sion in JHK bands of 6611-20515 not coming from the disk is 
that its QjHHK index, which is defined without using any opti- 
cal band, is entirely compatible with a reddened photosphere. 
Finally, in the left panel of Fig.|2]it is evident that the expected 
photospheric flux is larger than the observed one, since the disk 
is so inclined to partially occult the central star. We conclude that 
this star, and the other 1 1 sharing the same characteristics, does 
not show any NIR excesses in the IRAC and bands, but instead 
a V excess induced by scattered light. 



The first sample is composed of 90 stars with excesses in in- 
frared colors and with V - I bluer than the cluster PMS locus in 
the y vs. y - / optical diagram. Thus, they are apparently older 
than the adopted upper limit of the cluster age. We discuss pho- 
tometric and spectroscopic evidence supporting the membership 
of these stars to NGC 661 1. SED analysis suggests that their op- 
tical colors are altered by a significant amount of light scattered 
by circumstellar material into the line of sight, which rises the 
observed optical flux in V band. We present hints that this effect 
is more evident in young stars (younger than lO*" years) prob- 
ably as a result of a higher density of small grains in the disk 
atmosphere than in disks around older stars. In several cases, a 
partial occultation of the stellar photosphere by an inclined disk, 
alone or together with scattering effects, can be invoked. The oc- 
cultation, in fact, decreases the flux observed from the sources, 
but the optical colors remain unchanged. We conclude that the 
effects related to the presence of a disk (mostly if it is observed 
at high inclination) shift the position of disk-bearing stars in the 
y vs. y - / diagram from the PMS to the Main Sequence locus, 
altering the estimate of stellar mass, effective temperature, and 
age. 

The second sample of stars is composed of 66 stars with NIR 
excesses and photospheric IRAC colors. Our analysis suggests 
that these stars are disk-bearing members with 

- a disk with a large inner hole and a broad siUcate band in 
emission: in these cases the emission from the disk becomes 
significant only at 8.0//m and at longer wavelengths; 

- or a disk with a low mass (i.e. with MdisklMsmr < 10"^) pro- 
ducing small excesses in NIR bands. These excesses are not 
evident in color-magnitude diagrams, in which these stars 
can be easily confused with normal colors stars, but they are 
selected by our peculiar disk diagnostic, which compares in- 
frared and optical colors; 

- or a disk that significantly scatters the stellar optical light 
into the line of sight. These stars are not selected thanks to 
excesses in infrared bands induced by the disk, but thanks 
to the light scattered by the disk, which affects the optical 
colors. 
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6. Summary 

In this paper we examine the SEDs of two samples of pe- 
cuhar stars with excesses in 2MASS, UKIDSS, and IRAC 
bands, which were selected in our previous papers devoted to 
NGC 661 1 and M16. The excesses were detected using suitable 
reddening-free color indices. 
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Fig. 7. The observed SED (points) of the star 6611-20515. In the left panel the 5 compatible models (lines) and the predicted 
photospheric flux (dashed line) are also shown. The right panel shows the components of the best-fit model. As in Fig. [3] all the 
SEDs are normalized at the distance of 1 Kpc, and the observed SED in the right panel has been dereddened using the visual 
extinction estimated by the SED model (Ay = 2.24"'). 
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